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ABSTRACT 

Several classes of neutron star are sources of coherent emission at frequencies of 10^ — 
10^ MHz: others are radio-quiet. The primary emission spectra are broadly universal 
in form over many orders of magnitude in rotation period and polar-cap magnetic 
flux density. The existence of nulls and mode-changes in some radio-loud pulsars can 
be understood only as a manifestation of magnetospheric bistability. An ion-proton 
plasma with a possible background of electron-positron pairs is formed at the polar 
caps of stars with positive corotational charge density and is shown here to be a 
physical basis for the presence or absence of coherent emission and a likely reason why 
bistability may be present in the later stages of a pulsar lifetime. 
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1 INTRODUCTION 

Coherent emission by pulsars usually has a flux density 
Slv) oc in the frequency interval 10^ < u < 10® MHz 
with a typical spectral index a ~ 1.8. It is remarkable that 
this is true over intervals of three orders of magnitude in the 
rotation period P and five orders of magnitude in the spin- 
down inferred magnetic dipole field Bd- It is also remarkable 
that the luminosities estimated by Szary et al (2014) from 
the flux densities 5'i4oo for the complete ATNF pulsar cata¬ 
logue (Manchester et al 2005) are uncorrelated with position 
in the P — P plane. It is now widely accepted that the source 
is a plasma instability, possibly plasma turbulence, at low 
altitudes above the polar caps within the open sector of the 
magnetosphere. Exceptional instances have been observed in 
which the spectral index is a ~ 0 or in which the emission is 
both intermittent and outside the normal well-defined win¬ 
dow of pulse longitude. We shall refer to these as secondary 
coherent emission and it is possible that new high-sensitivity 
systems, such as the Square Kilometre Array (SKA) will lead 
to the discovery of further examples. 

Coherent emission at radio frequencies requires either 
maser action between suitable energy levels or macroscopic 
fluctuations in charge and current density. Unlike incoher¬ 
ent processes such as the emission of X-rays or y-rays, the 
nature of the charged particles contained in these fluctua¬ 
tions is not of immediate significance for the actual radiation 
process. Emission can depend only on the form of the fluc¬ 
tuations. The leading incoherent processes must certainly 
be those of lowest-order quantum electrodynamics involving 
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electrons or positrons. But in the case of coherent emission, 
there is no requirement that the particles involved should be 
electrons or positrons. Sources of incoherent emission are, 
with few exceptions, thought to be at high altitudes in the 
magnetosphere, near the light cylinder, but are not the sub¬ 
ject of this paper. 

The existence or otherwise of a suitable plasma insta¬ 
bility does, however, depend on the plasma composition and 
we do not believe that progress can be made without due 
attention to this. The relation between rotational angular 
velocity fl and polar-cap magnetic flux density B is of im¬ 
mediate significance. For the case of fl ■ B >0, the polar-cap 
Goldreich-Julian charge density is negative. The accelerated 
particles must be electrons and electron-positron pair cre¬ 
ation the only possible secondary process: a reverse flux of 
positrons has no effect apart from heating the polar cap. 

The Q-B < 0 case is much more complex. A reverse flux 
of electrons creates protons which are an important compo¬ 
nent of the accelerated plasma. Here, we feel obliged to say 
that empirical methods are not apt for the determination 
of plasma composition. The nuclear physics processes are 
well-established and given their acceptance, there can be no 
serious doubt that a proton component is present. As a con¬ 
sequence of the form of the electromagnetic cascade shower 
produced by reverse electrons, two time constants are intro¬ 
duced, both with observable consequences. 

Plasma properties in the f2 ■ B < 0 case have been 
investigated in a number of previous papers (Jones 2010; 
see also Jones 2013a and papers cited therein). Further pa¬ 
pers (Jones 2013b, 2014a,b) based on this work gave rea¬ 
sons for believing that the secondary electron-positron pair 
plasma which may be formed in pulsars with either sign 
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of 17 ■ B is not the source of the observed primary coher¬ 
ent emission. Given the ineluctable evidence that the basic 
spectral properties of the primary coherent radiation are 
universal, the reluctance to consider sources other than a 
secondary electron-positron plasma is surprising, although 
this system might be relevant to processes such as reso¬ 
nant cyclotron emission (Kazbegi, Machabeli & Melikidze 
1991) or cyclotron-Cherenkov resonance (Lyutikov, Bland- 
ford & Machabeli (1999) occurring near the light cylinder. 
But neither these sources of light-cylinder coherent emission 
nor other electron-positron sources such as collective curva¬ 
ture radiation (Istomin, Philippov& Beskin 2012) and the 
Langmuir instability in a non-stationary plasma (Asseo & 
Melikidze 1998) will be considered further in this paper. 

The previous series of papers referred to above makes no 
specific assumption about the nature of the global magne¬ 
tosphere and in this respect is at variance with much recent 
work which uses numerical plasma-physics techniques (for a 
comprehensive list of published papers see Tchekhovskoy, 
Spitkovsky & Li 2013). Model force-free magnetospheres 
have been constructed extending outwards from a typical 
inner radius of 0.2Rlc with a boundary condition allowing 
particles of either sign to pass freely through to regions be¬ 
yond the light cylinder radius Rlc- In this work, the global 
magnetosphere (outside the inner radius of the model) de¬ 
termines the current density in the open magnetosphere. 
The charge density differs little from the Goldreich-Julian 
value, but the current density can be significantly different, 
either larger or of the opposite sign. Such a state can be 
formed only by the counterflow of oppositely charged parti¬ 
cles which, in the $7 • B > 0 case usually considered, must 
rely on electron-positron pair creation. The consequences of 
a constrained current density for the form of magnetosphere 
at low altitudes above the polar cap have been investigated 
recently by Ghen & Beloborodov (2013), Timokhin & Arons 
(2013) and Philippov &Spitkovsky (2014). Pair creation has 
been modelled in detail by Timokhin & Arons who adopted 
the not unusual assumption of an extremely small ( 10 ® cm) 
flux-line radius of curvature. The present paper is not a cri¬ 
tique of this work; it makes the different assumption that 
the current density is determined at low altitudes above 
the polar cap. It emphasizes the importance, in the first 
place, of establishing the physics of the polar-cap region of 
high energy and particle number densities under conditions 
such that counterflow at fluxes of the order of the Goldreich- 
Julian density is not possible, in which the magnetic field B 
is almost precisely immutable, and particle motion is parallel 
or anti-parallel with it. 

The form of global magnetosphere which would exist 
given the flux of relativistic particles in the open sector, 
with moderate Lorentz factors of the orders of 10^“^, as a 
boundary condition at 0.2Rlc is not known. Observational 
tests of its nature are not obvious in a normal (middle-aged) 
pulsar. Goherent emission is from the polar cap region, levels 
of incoherent emission from the vicinity of the light cylin¬ 
der are not usually high, and the spin-down index probably 
unobservable owing to noise. In any case, the connection be¬ 
tween any of these observables and global magnetospheric 
structure is by no means direct. 

Our proposal here is that the primary coherent emis¬ 
sion has a common origin (as its almost universal form ap¬ 
pears to indicate) in multi-component Langmuir modes in 


the magnetospheres of ^7 • B < 0 neutron stars. This was 
discussed in Jones (2014b) but that paper had the serious 
deficiency of failing to describe adequately how the transi¬ 
tions occur in what can be a bistable magnetosphere with 
nulls and mode changes. It has since been found (Jones 
2014c) that the dielectric tensor for an ion-proton plasma 
which also includes a background component of electrons 
and positrons gives a zero Langmuir-mode growth rate for a 
certain class of electron-positron number densities and spec¬ 
tra. This can be the case even for densities small compared 
with the Goldreich-Julian value, and is no more than a con¬ 
sequence of the small electron-proton mass ratio. It appears 
that this property of multi-component modes allows a mod¬ 
est change in the electron-positron spectrum to extinguish a 
Langmuir mode which would otherwise have a high growth 
rate. We shall attempt to describe why such changes can 
occur naturally, with certain time-scales, in Section 2.2. 

A more extensive summary of the plasma properties 
is contained in Sections 2.1 and 2.2. The polar-cap system 
described there is complex with many variables in compari¬ 
son with the opposite case, 17 ■ B > 0, for which Maxwell’s 
equations, the processes of lowest-order quantum electro¬ 
dynamics, and the laws of mechanics suffice. Its properties 
also depend, unfortunately, on two parameters fixed for a 
given neutron star which are not well known, specifically, 
the whole-surface temperature Ts and the atomic number 
Z® of nuclei in the thin surface layer which is the source 
of ions at the polar cap. As a consequence of the complex¬ 
ity, the system in its general form has not been completely 
modelled, but the qualitative results of the limited modelling 
that has been done are described. 

This paper examines a limited question. Are the proper¬ 
ties of the ^7•B < 0 plasma which are summarized consistent 
with the broad properties of the coherent emission and its 
presence or absence in the various classes of pulsar or iso¬ 
lated neutron star that are observed? It is not concerned 
with the sui generis cases that are not uncommon in pulsar 
observations or with more detailed considerations such as 
the interpretation of polarization which may be intrinsic or 
the result of propagation through the magnetosphere. The 
effects of the latter may be of significance for estimates of 
the radiation emission height but are not considered here. 
The question of consistency is addressed in Section 3, with 
some brief notes on the problems that remain being given 
in Section 3.6. Our conclusions are stated in Section 4. It is 
hoped that clarification of polar-cap physics will be useful in 
relation to work on the incoherent processes near the light 
cylinder which are the source of X-ray and 7 -ray emission. 


2 PHYSICS OF THE POLAR CAP 

A qualitative understanding of how plasma composition de¬ 
pends on neutron-star parameters is necessary for the pur¬ 
poses described in Section 1. In order to achieve this, a 
fairly comprehensive summary of some of the conclusions 
reached in previous papers (Jones 2013a and earlier papers 
cited therein) with a statement of the physical factors in¬ 
volved is required and is the subject of this Section. A clear 
distinction between negative and positive Goldreich-Julian 
polar-cap charge densities is necessary. Our summary refers 
to positive densities unless otherwise stated. 
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2.1 Plasma composition 

The basis for quantitative calculations has been a circular 
polar cap of radius uo satisfying space-charge limited flow 
boundary conditions. Thus there is an atmosphere in local 
thermodynamic equilibrium (LTE) of which electrons and 
ions are the dominant components. Its scale height is of the 
order of 10~^ cm, many orders of magnitude smaller than 
Mo- A density discontinuity exists between the atmosphere 
and the solid (or liquid) phase of the surface. The ion density 
at the foot of the atmosphere is not well known but is likely 
to be of the order of 10^^ — 10^"^ cm“® . 

A reverse flux of electrons from pair creation or photo¬ 
ionization must be incident on the polar cap. The electro¬ 
magnetic showers formed in the atmosphere, or in the at¬ 
mosphere and solid surface, produce protons at rates which 
have been calculated and are linear functions of shower en¬ 
ergy. The electric field within the LTE atmosphere is deter¬ 
mined principally by the electrons and ions and is such that 
the relatively small density of protons cannot be in static 
equilibrium: they experience an upward force to the top of 
the atmosphere. This fractionation is a direct consequence 
of the difference in charge-to-mass ratio between ions and 
protons. Therefore, the Goldreich-Julian flux is formed pref¬ 
erentially from protons, surplus protons forming a separate 
atmosphere above the ions. It has an ion component only at 
times when the proton atmosphere has been exhausted. The 
diffusion time, from the point of formation in the shower to 
the top of the LTE atmosphere has been estimated to be of 
the order of Tp ^ 1 s (see Jones 2013a). 

During the lifetime of a pulsar, an extremely thin layer 
of mass no more than about 10 ~® Mq over the whole surface 
of the star is sufficient to supply an ion flux at the Goldreich- 
Julian density. This is so small that its atomic number 
is likely to depend on unknown details of the formation pro¬ 
cess such as fallback and so must be regarded as an un¬ 
known parameter. The lateral diffusion coefficient for pro¬ 
tons, whether derived from the Goulomb collision frequency 
or from the Bohm mechanism (see Montgomery, Liu & Va- 
hala 1972) is small, no more than of the order of 10“^ —10~^ 
cm^ s“^, so that the outward proton flux at any point u on 
the polar cap is determined only by the inward electron en¬ 
ergy flux at previous times at that point. 

Proton (and neutron) production in an electromagnetic 
shower is almost entirely through the formation and decay 
of the nuclear giant-dipole state by photons in the 15 — 30 
MeV energy interval. Consequently, it is localized at the 
shower maximum which is typically at about 10 radiation 
lengths depth. The properties of a shower at B ~ 10^^ G are 
not essentially different from the zero-field case, although at 
fields greater than the critical Be = 4.41 x 10^® G there 
are significant differences, among them being the increased 
importance of the Landau-Pomeranchuk-Migdal effect. Pro¬ 
ton production rates have been estimated in this case but 
the differences are not large enough to change any of this 
paper’s conclusions.. 

Photo-disintegration through the giant-dipole state re¬ 
duces the nuclear charge from Z^ to Zs as ions move up¬ 
wards to the polar-cap surface. There is a delay between 
proton formation and its effect on the nuclear charge of ions 
accelerated from the surface. This introduces a second time 
constant of interest, specifically, the time Tri in which ions 


equivalent to one radiation length leave the surface. Both Tp 
and Tri are consequences of shower structure and both lead 
naturally to instabilities in the composition of the acceler¬ 
ated plasma. 

Ions of nuclear charge Z^ have an ionic charge Z in the 
LTE atmosphere which is dependent on the polar-cap tem¬ 
perature Tpe- There is a small degree of further ionization 
Z ^ Zh ss ions move initially through the inertial acceler¬ 
ation region (Michel 1974) of the polar-cap blackbody ra¬ 
diation at very small altitudes z uo above the polar-cap 
where ion Lorentz factors are no more than of the order of 
unity. Further acceleration is a consequence of the Lense- 
Thirring effect (Muslimov & Tsygan 1992). But the more 
significant ionization is at higher altitudes and Lorentz fac¬ 
tors in the acceleration region 1 < rj < rja, where rja is 
a suitable upper limit to the acceleration region. Photons 
from a large area of the surface at temperature Ts can be 
boosted in this region by Lorentz transformation to energies 
sufficient for the ionization of more deeply bound electrons. 
(Here the polar coordinate is expressed as r = rfR in terms 
of the neutron-star radius R.) The final ion charge is Z^o- 
Tables of Z^o and of the mean reverse-electron energy per 
ion e have been calculated as functions of the acceleration 
potential difference Voo for that ion, treated as an indepen¬ 
dent variable. These have been obtained for various B, Za, 
and for whole surface temperatures Ts = 1 — 4 x 10® K . 

In terms of electrostatics in the corotating frame of 
reference, photo-electrons accelerated inward to the surface 
have precisely the same effect as pair creation in partially 
cancelling the acceleration field. The electron charge density 
distribution required to cancel the acceleration field exactly 
can be found directly from the Lense-Thirring effect expres¬ 
sion for the Goldreich-Julian charge density aojitf) (Mus¬ 
limov & Tsygan 1992, see also Harding & Muslimov 2001). 
Its density is. 


(Tein) = 


K rf 

2rp 




( 1 ) 


and the mean electron track length is R/2. Here, k. ~ 0.15 is 
the Lense-Thirring factor which is a function of the neutron- 
star radius and moment of inertia, but is independent of B 
and P. Polar-cap modelling (Jones 2013a) has shown that 
the maximum acceleration field which would be present in 
the absence of photo-ionization or pair creation is largely 
cancelled by the reverse electron flux. Thus It is a fair ap¬ 
proximation to assume that the electron track-length distri¬ 
bution is a universal constant. Photo-ionization transition 
rates are dependent on the ratio of the threshold photon 
energy in the ion rest frame to the blackbody photon en¬ 
ergy at temperature Ts. Hence it follows that Ion Lorentz 
factors and reverse-electron energies are oc approxi¬ 

mately. This is true for the mean values of these quantities 
even though the modelling shows that the electric field and 
plasma composition fluctuate over the polar-cap. This result 
is of significance in understanding why coherent emission 
ceases late in the lifetime of a pulsar. 

If the magnetic flux density is small, as in the case 
of millisecond pulsars (MSP), the process is effectively un¬ 
changed because zero-field photo-electric cross-sections are 
adequate. As we have already mentioned, shower proper¬ 
ties are also essentially unchanged. The role of the magnetic 
flux density in the process is merely to provide an acceler- 
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ation potential above the polar cap of a rotating neutron 
star. We emphasize that, in this respect, the open magneto¬ 
sphere described here differs completely from one in which 
electron-positron pair creation is the essential process and, 
in particular, from those of f2 ■ B > 0 neutron stars. 

The state of the polar cap is chaotic as a consequence 
of the local nature of the shower processes. This has been 
confirmed by the construction of a self-consistent discrete¬ 
cell polar-cap model (Jones 2013a) based on the tabulated 
data described above which gave evidence of conal emis¬ 
sion and, though limited by the assumption of constant Zs, 
demonstrated fluctuations in plasma properties not inconsis¬ 
tent with the subpulse modulation that is observed almost 
universally in the coherent radio emission (Weltevrede, Ed¬ 
wards & Stappers 2006). But even though Zg is likely to be 
a constant over the polar-cap area, there is no reason to as¬ 
sume that the atomic number of nuclei reaching the surface 
will be other than some function Zs{u,t) of no simple form 
in either position or time. In this respect, our modelling has 
been incomplete. 

2.2 Langmuir modes 

The physics considerations summarized above are individ¬ 
ually well understood and we regard them as establishing 
the nature of the accelerated plasma in ■ B < 0 neutron 
stars beyond any reasonable doubt. The generation of co¬ 
herent emission is extremely efficient (see Jones 2014b) and 
the well-defined velocity distribution of an ion-proton sys¬ 
tem indicates that Langmuir modes should be examined as 
a possible source. The longitudinal mode does not couple 
with the radiation field, but Asseo, Pelletier & Sol (1990) 
demonstrated that quasi-longitudinal modes can have simi¬ 
lar growth rates and do couple effectively with modes that 
are capable of leaving the magnetosphere. However, there 
has since been interest in the generation of plasma tur¬ 
bulence from the longitudinal mode (see Weatherall 1997) 
and it is not unreasonable to suppose that this is possible 
through Langmuir-mode growth in a plasma which is not ho¬ 
mogeneous and is laterally of finite dimensions. Under these 
conditions, we shall assume that coupling with the radiation 
field does occur, converting particle longitudinal kinetic en¬ 
ergy to radio-frequency emission. 

Langmuir-mode angular frequencies for the ion-proton 
plasma tend to be no more than of the order of the low- 
frequency turnover that is observed in some pulsar spectra. 
This is not a problem if we assume (as in Jones 2013b) that 
the development of turbulence transfers energy to higher 
wavenumbers as is the case, for example, in homogeneous 
isotropic fluid turbulence. However, no quantitative studies 
appear to have been made for the present case, although 
the numerical plasma methods of Weatherall (1997, 1998) 
have shown that, in electron-positron plasmas, very high 
wavenumbers can be generated such as are observed in the 
Crab pulsar (Hankins & Eilek 2007). 

The maximum particle longitudinal kinetic energy that 
can be transferred by any mechanism to electromagnetic ra¬ 
diation is an important quantity. Luminosities estimated for 
a group of normal (middle-aged) pulsars by Jones (2014b) 
showed that the radio-frequency energies generated are typi¬ 
cally in the interval W — 1 — 10 Gev per unit positive charge 
in the primary beam, assuming a Goldreich-Julian density. 


For the ion-proton beam considered here, the kinematic- 
predicted value of W is given by equation (14) of that paper, 

{1 + y){x'yp + yji) 

where x = A/Z is the mass-to-charge ratio of the ion, y 
is the mean number of free protons in the beam per unit 
positive charge of the ions, are the Lorentz factors and 
nip is the proton mass. Its dependence on baryon masses 
indicates that W can be of the observed order. 

Many of the views expressed in this present paper are 
based on recent studies of Langmuir-mode growth rates for 
ion-proton plasmas which have an additional small electron- 
positron component (Jones 2014c). There are several possi¬ 
ble sources of pairs other than the curvature radiation (CR) 
assumed in the classic model of Ruderman & Sutherland 
(1975) and in many later papers. The most quoted mecha¬ 
nism is inverse Gompton scattering (IGS) in which a black- 
body photon is scattered by an outward moving electron, 
with magnetic conversion of the resultant photon of high 
momentum k when its transverse component exceeds the 
threshold k±_ > 2mc. (Directions are with respect to the 
local B.) The polarization-averaged photon attenuation co¬ 
efficient in the classical limit defined by k±_ ^ 2 mc and 
X = k±B/2mcBc <C 1 is, 

^ exp(-4/3x)) , (2) 

(Erber 1966). Gonversion lengths ~ uo are given by x 0.1, 
which is in the classical region and can be attained for 
B Be- The IGS process can be significant for 12 • B > 0 
neutron stars (see Hibschman & Arons 2001, Harding & 
Muslimov 2002) but for the opposite sign, it must be con¬ 
siderably less effective because there is no primary source of 
outward-moving positrons, with the consequence that any 
high-energy Gompton-scattered photons enter the polar-cap 
surface without gaining the pair-threshold momentum per¬ 
pendicular to B. Electromagnetic showers also produce neu¬ 
trons, a small fraction of which may diffuse back to the top 
of the LTE atmosphere. Therefore, photons from (n, 7 ) re¬ 
actions near the top of the atmosphere are a further source 
but its importance is difficult to estimate quantitatively. In 
common with IGS, this process requires a polar-cap field 
strong enough for single-photon magnetic conversion. The 
only source not essentially dependent on B appears to be 
pair production by photon-photon scattering, but this can 
be significant only for very high polar-cap temperatures. 

The most important result of the study was that the 
presence of an electron-positron component sparse compared 
with the Goldreich-Julian number density can reduce the 
longitudinal Langmuir-mode growth rate to zero if its ve¬ 
locity spectrum overlaps the ion and proton velocities. Gon- 
sequently, a small change in the characteristics of a sparse 
electron-positron component can effectively switch off the 
turbulence. An interesting question which then arises is 
whether or not the whole polar cap is likely to be affected. 
It has been shown (see Fig. 1 of Jones 2014c,) that, for 
a uniform plasma characterized by a given set of parame¬ 
ters, there is a wavenumber of maximum growth rate. But 
if we regard the wavenumber as an independent variable, 
the mode has growth rates of a similar order over an inter¬ 
val that is finite but not unlimited. It is undoubtedly the 
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case that the polar-cap plasma is of finite lateral dimen¬ 
sions and is far from homogeneous in either composition or 
Lorentz factor. (We assume a zero-potential boundary con¬ 
dition on the cylinder of radius uo{rj) separating open from 
closed magnetosphere.) 

Under these conditions, what is the state of the po¬ 
lar cap at some general instant of time? The general free- 
vibrational state of an inhomogeneous elastic body would 
simply be a superposition of its normal modes, though ob¬ 
taining these does present problems (see, for example, Viss- 
cher et al 1991). This suggests that, in the present case, we 
should consider the growth of a Langmuir mode from an 
initial fluctuation that is confined to a localized area of the 
polar cap. On the basis that growth rates exist over a finite 
interval of wavenumber in a homogeneous plasma, it is likely 
that a localized mode would exist with significant amplitude 
gain on flux lines passing through that area even though it 
encompasses some variation in plasma parameters. The gen¬ 
eral state of the polar cap would then be a superposition of 
such modes. If this is accepted, we can see that a modest 
change in the electron-positron component in some area of 
the polar cap can effectively switch off a mode. In this, we 
can see the basis for both nulls and mode-changes which will 
be described in the following Section. 


3 A BASIS FOR UNDERSTANDING THE 
OBSERVED COHERENT EMISSION 


The polar-cap physics of D ■ B <0 neutron stars is necessar¬ 
ily dependent on parameters that are not well known. These 
are, principally, the polar-cap surface atomic number and 
magnetic flux density B, and the whole-surface blackbody 
temperature Ts . It has been assumed here and for modelling 
purposes that B can be approximated by the ATNF pulsar 
catalogue values Bd, there being no alternative (Manchester 
et al 2005: in modelling the polar cap we have assumed a 
radius R = 1.2 x 10® cm, which is larger than the ATNF 
catalogue assumption). The whole-surface temperature we 
refer to here is, more precisely, that for a circular area cen¬ 
tred on the polar-cap, subtending of the order of a steradian, 
but with photons in a frame intermediate between the ob¬ 
server and the proper frame of the neutron-star surface. The 
most recent survey of existing neutron-star surface temper¬ 
atures has been published by Keane et al (2013) who also 
add three upper limits. Temperatures listed for ages up to 
10 Myr are not inconsistent with the interval of T, used in 
modelling the polar-cap (Jones 2013a). There appears to be 
no direct evidence on the third parameter, Z®, which is not 
unconnected with T,. Ions with very low Za have relatively 
weakly-bound electrons and hence charges Z — Zoo = Zs- 
In this case, there would be no photo-ionization and there¬ 
fore no reduction in the acceleration field and in ion Lorentz 
factors, which in most neutron stars would reduce the mode 
growth-rate to insignificant values. The mode amplitude is 
oc expA(r;), in which. 


A = 


RT 


c 



* T / 

uj^drj 


(3) 


approximately, at altitude rj, where T is the dimensionless 
growth-rate given in Fig. 1 of Jones (2014c) whose value is a 
slowly-varying function of the relative ion and proton fluxes. 


and uj* is defined by. 


‘iTvNo{ri)Z^e^ 


(4) 


in terms of the mass Mz, number density Nz and Lorentz 
factor 7 (r;) of the accelerated ion. In general, the most 
significant contributions to A are from small rj, where 7 
is also small. Assuming a constant 7 , the growth-rate de¬ 
pendence on the basic neutron-star parameters is given by 
A oc {—B excluding the parameter T 

and the charge-to-mass ratio of the ion in the LTE atmo¬ 
sphere. (Here ip is the angle between B and fl). It has 
been assumed in previous papers that A > 30 is adequate 
for the development of turbulence. Thus even if a proton 
component were present, acceleration of very small-Zg ions 
to high Lorentz factors would not allow adequate growth 
of the Langmuir mode. The processes described in Section 
2.2 would not function. But provided this limit is not ap¬ 
proached, the value of Z® is not critical. 


3.1 Factors essential for coherent emission 

With these qualifications, we can see qualitatively how the 
presence or absence of coherent emission depends on B and 
Ts. In comparisons with observation, we can assume only 
that Ts decreases as a function of characteristic age U = 
P/2P. Then the polar-cap modelling in Jones (2013a) and 
the Langmuir-mode studies (Jones 2014c) reveal two distinct 
factors. 

(i) Photo-ionization transitions require black-body pho¬ 
ton momenta khb to exceed the threshold, immediately above 
which cross-sections are large. Very approximately, they 
should satisfy 'ykhbC = Eb, where Eb is the separation en¬ 
ergy of the least bound electron at ion charge Z. Thus the 
Lorentz factor required for photo-ionization is 7 oc TpcTp^. 
Hence decreasing Ts at constant B requires that ions be 
accelerated to higher Lorentz factors for photo-ionization, 
as shown in Section 2.1. As the limit T, ^ 0 is approached, 
the reverse-electron energy per ion (e) first increases, oc Tp^, 
as does the rate of proton production. But if ions have to 
be accelerated to very high Lorentz factors for the required 
photo-ionization rate, the mode-amplitude growth exponent 
A given by equation (3) will be decreased to the extent that 
growth-rates in the limit fall to values too small to allow the 
development of turbulence and observable coherent emis¬ 
sion. The pulsar then enters a permanent radio-quiet state. 

(ii) The rate of ICS pair creation depends on both B 
and Ts- There are two separate reasons. Equation (2) shows 
that single-photon magnetic conversion at 10^^ G requires 
a transverse photon momentum k± « 9mc. But at much 
higher fields B ~ Bo, pair creation transition rates are high 
immediately above the threshold fcx = 2mc. Whilst trans¬ 
verse momentum is not conserved at finite B, the final ICS 
photon fcx must be balanced by a corresponding component 
in the Fourier transform of the electron Landau function. 
Thus the probability that an ICS black-body photon would 
reach the threshold for pair creation is small at 10^^ G but 
increases as B ^ Bo- Hence the probability of mode sup¬ 
pression increases rapidly as a function of B. Secondly, the 
kinematics of ICS in the Klein-Nishina region of energy for 
an inward-moving electron show that the threshold energy 
for production of a final-state photon with k± > 2 mc is 
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(X T~^. the blackbody photon density is oc T® so that the 
probability of pair creation at a fixed is a rapidly increas¬ 
ing function of Ta. An increased electron-positron density 
with the appropriate momentum spectrum tends to suppress 
one (or more) Langmuir modes and hence the turbulence de¬ 
veloping on a given set of flux lines according to arguments 
given in Section 2.2. 

Factors (i) and (ii) determine the conditions under 
which coherent emission is possible. But in using them we 
have to remember that B may differ significantly from Bd 
and that Zg is unknown. 

Some of the relevant observational data concerning dif¬ 
ferent sets of pulsars are shown in Fig.l. These are not in¬ 
tended to be complete but are representative of the class 
concerned. On the basis of equation (2) and of the char¬ 
acteristics of curvature radiation, the pair-creation thresh¬ 
old can be approximated by a hxed value of the variable 
X = Bdi 2 P~'^^'^, assuming a dipole field (see also Fig. 1 of 
Harding & Muslimov 2002). A threshold value Xc = 6.5 is 
consistent with Harding & Muslimov. 

Normal pulsars are too numerous to be shown individu¬ 
ally: the vertical broken lines show the half-maxima of their 
distribution in Bd and most lie at A < Xc- Pulsars exhibit¬ 
ing nulls and nulls plus mode changes are shown separately, 
and are those listed in Tables 1 and 2 of Wang, Manchester 
& Johnston (2007). Radio-loud millisecond pulsars (MSP) 
are also shown in separate groups: those listed in the second 
Fermi LAT catalogue (Abdo et al 2013) and those which 
have no observable y-emission but are contained in Table 
3 of Kramer et al (1999). The Rotating Radio Transients 
(RRAT) are those listed with timing solutions by Keane et 
al (2011). The Central Compact Objects (CCO) are those 
described by Gotthelf et al (2013): the X-ray Isolated Neu¬ 
tron Stars (XINS) and High-Field Rot at ion-Powered Pulsars 
(HBRPP) are those listed by Kaplan & van Kerkwijk (2009). 
The magnetars are from the McGill catalogue (Olausen & 
Kaspi 2014). 

3.2 Normal pulsars 

Pulsars not placed in the specific classes listed here are re¬ 
ferred to as normal: they are radio-loud and comprise the 
major part of the ATNF catalogue. The typical field is 
B ~ 10^^ G so that the classical limit for single-photon 
magnetic conversion is apt. A small fraction of them lie at 
X > Xc and the immediate question is whether in these 
there is self-sustaining CR electron-positron pair creation at 
the polar cap. Any having ■ B > 0 should produce pairs 
but our view (Jones 2013b, 2014b) is that they would not 
be radio-loud in the sense of having the usual large negative 
spectral index. It is possible to argue that a state of pair 
formation at A > Ac in those with f2 • B < 0 would be un¬ 
stable against a transition to the ion-proton state. Ion emis¬ 
sion with formation of reverse photo-electrons must always 
be possible in polar-cap areas near the boundary uo where 
the acceleration field is not large enough to support CR pair 
creation. Any inward fluctuation of this zone must further 
reduce the acceleration field and so enlarge its area. For this 
reason, we believe that an ion-proton plasma is the source 
of coherent emission even in very young pulsars. Polar-cap 
modelling (Jones 2013a) has shown, unsurprisingly, that rel¬ 
atively high surface temperatures, Ts ^ 5 x 10® K, much re¬ 


duce both the mean acceleration potential difference above 
the polar cap and fluctuations away from this state, which 
is therefore effectively stable. 

Consideration of factors (i) and (ii) for normal neu¬ 
tron stars at A < Ac with $7 • B < 0 and Bd 10^^ G 
shows that ICS pair production is negligible compared with 
the Cl ■ B > 0 case. Also the polar-cap temperature Tpc is 
not high enough to give significant rates for pair produc¬ 
tion by photon-photon scattering. Thus for the ■ B < 0 
case, with moderate values of J3, a radio-loud state should 
be expected until the final turn-off. Observable coherent 
emission ends because the amplitude growth rate exponent, 
A oc B^^^, derived from equations (3) and 
(4), must eventually become small as Ta declines. Here, we 
neglect cos ip, also any -i/i-dependence in our use of the ATNF 
catalogue values of Bd- 

The heat flux represented by Tpc has two sources: the 
reverse-electron energy and the heat loss from the interior of 
the star which is the source of Ta - The time-averaged reverse- 
electron energy flux has an upper limit equal to that which 
would create a Goldreich-Julian flux of protons (see Jones 
2010; equation 33). The shower energy required to produce 
one proton is a slowly-varying function of B but is treated 
here as a constant. For high surface temperatures exceeding 
^ 100 eV, it is easy to confirm that for a normal star the heat 
flow from the interior exceeds that of the reverse flux. Then 
Tpc = Ts is a fair approximation. But the eventual decline 
of Ta and the increase in P will reduce A so that growth 
becomes too small for the development of turbulence. 

The presence of B in equation (3) derives from the 
Goldreich-Julian charge density: it is therefore the true mag¬ 
netic flux density, not Bd- In principle, B — (pBd with C, 
larger or smaller than unity: an example would be the field 
of a non-central dipole. Therefore in Fig. 1, which gives 
as a function of Bd for the various groups con¬ 
sidered here, it should not be expected that a constant B- 
independent cut-off value will be seen. 

In fact, it is possible to superimpose on Fig.l a curve 
showing a well-defined cut-off over the first five orders of 
magnitude in Bd- The multitude of normal pulsars in the 
ATNF catalogue is too large to be shown, but the relevant 
boundary of their distribution has a sharp cut-off shown in 
the diagram as an arc of a curve. It would be consistent with 
the imaginary superimposed curve referred to above, were 
we to draw it. The expression for A does not necessarily 
apply for Bd > Be because mode growth-rates are liable to 
be reduced to zero by pair creation, as described in Section 
3.1 (ii). 

The indicated cut-off value of By^^P~^^^ is not precisely 
independent of Bd but varies remarkably slowly over the five 
orders of magnitude to which we refer. The expression for 
A has been obtained by the simplest scaling considerations 
and too much should not be expected of it. It is also possible 
that there are systematic rather than random variations in 
(P- Thus (( > 1 at small Bd (the MSP and GGO) would render 
the true cut-off less B-dependent. It may be relevant that 
most MSP are in binaries, and their history is quite different 
from that of the normal multitude. 

We have noted previously, in Section 2.1, that the fluc¬ 
tuations in polar-cap acceleration potential and plasma com¬ 
position that are the result of the time constant Tp are not in¬ 
consistent with the sub-pulse modulation which is observed 
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Figure 1. The factor shown as a function of for the sets of pulsars and neutron stars listed in the final paragraph 

of Section 3.1. These are not in all cases complete lists of their class but are convenient representative groups. The large number of 
normal pulsars in the ATNF catalogue are not shown individually but vertical broken lines show the half-maxima of their distribution 
in From the ATNF catalogue, it has been possible to deduce the lower boundary of their distribution and this is shown as a solid 
arc between the half-maximum limits. The MSP are divided into those with a Fermi LAT catalogue entry (solid circles), all of which 
are radio-loud, and a radio-loud set with no observed 7 -emission. Pulsars showing nulls and nulls plus mode changes are also shown 
separately. 


almost universally. But it is likely that some of the modula¬ 
tion, particularly that observed at high time-resolution, is a 
consequence of the turbulent plasma state we believe to be 
the radiation source. 


Apart from the papers of Weatherall, which consid¬ 
ered pair plasma, there appears to be no work on radia¬ 
tive emission from the proton-ion turbulence we assume is 
formed. The nanosecond structure observed by Hankins & 
Eilek (2007) in the Crab pulsar main pulse provides possibly 
the most exacting test of any emission model. The example 
shown in Fig. 5 of their paper represents a nanopulse with 
an integrated flux of 2 x 10~^° erg cm“^. If the source at 
distance D is moving toward the observer with Lorentz fac¬ 
tor 7 s and, in its rest frame, emits into solid angle Si}s, the 
direction of emission in the observer frame can be at an 
angle of the order of 77 ^ with respect to B and the area 
illuminated is then ~ D'^SQs/Js ■ Given any reasonable as¬ 
sumptions about the two unknown parameters here, this is 
a remarkably large energy. Let As be the central wavelength 
emitted in the rest frame of the source. Collective particle 
motion and acceleration is parallel with B. Thus within the 
turbulent volume, a basic unit in the source rest-frame giving 
coherent emission might be a very short-lived quasi-dipole 
source, axis parallel with B and of length As, width wXs and 
depth \s/2'k. The energy available, per unit positive charge 
in the ion-proton plasma, for conversion to radio-frequencies, 
has been estimated previously (Jones 2014b) to be IT ~ 7 
Gev. The energy that can be extracted from this unit in 


principle is therefore no more than 
'wXl WBjy) 

27r7s Pec ’ ^ ’ 

in the observer frame. The coherent source volume can be 
large for w 3> 1, but the emission solid angle is then reduced 
to Jfis ~ 2-k/w. The insertion of parameter values for the 
Crab, with 7 s = 10, rj = 10, and the central frequency of 
9.25 GHz, gives As = 32 cm and w = 9x 10®. The total width 
w\s = 3 X 10^ cm is to be compared with an open magneto¬ 
sphere radius uo{ri) = rio(lO) = 2.8 x 10® cm. Clearly, some 
organization of coherence more subtle than our elementary 
example would be required. It must also satisfy the condi¬ 
tion that to produce a pulse length of 5t = 0.4 nanoseconds, 
the units must lie within a height interval of less than 2'yl cSt 
in the observer frame. Thus formation of nanosecond struc¬ 
tures is possible in principle, although we do not attempt 
here to estimate their probability. But it must be obvious 
how much more difficult it would be to understand these 
structures in terms of a light cylinder source. 


3.3 Millisecond pulsars 

Two sets of MSP are shown in Fig. 1; those listed in the 
Fermi LAT catalogue (Abdo et al 2013) and those not listed 
in that catalogue but present in Table 3 of Kramer et al 
(1999). A notable feature is the absence of radio-quiet MSP 
(defined as a flux density 51400 < 30 /r Jy) in the Fermi 
LAT catalogue. This is complicated by the different detec¬ 
tion procedures used; either known radio ephemerides or 
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blind searches using periodicities in LAT data. Nevertheless, 
Romani (2012) has estimated that the fraction of radio-quiet 
MSP can be no more than ~ 1/3. In view of the very small 
rotation periods and magnetic fields it is notable that the 
MSP radio spectra in both sets are very similar to those of 
normal pulsars (Jenet et al 1998, Kramer et al 1999, Es¬ 
pinoza et al 2013). 

A farther unresolved question concerns a small gronp of 
4 MSP in which the radio and 7 -emission light curves are in 
phase (Venter, Johnson & Harding 2012). The Crab pnlsar 
is the only other known example of this phase-alignment: 
phase leads or lags are otherwise the rule. Following the 
outer-gap model introduced by Cheng, Ho & Ruderman 
(1986) and further developed by Cheng, Ruderman & Zhang 
( 2000 ) to include pair creation by photon-photon collisions, 
it is now generally assumed that the source of 7 -emission 
must be close to the light cylinder. Phase-alignment has 
been used as evidence for coincident source regions, but this 
is open to question. 

With regard to the Crab pulsar, we have commented 
in Section 3.2 that it is not obvious how the nanosecond 
structure in its main pulse (Hankins & Eilek 2007) can be 
reconciled with emission from the vicinity of the light cylin¬ 
der. These indicate a polar-cap source for the coherent radio 
emission. If this is so, and the Crab is an $1 • B < 0 neutron 
star, the source of its generally accepted pair production 
must be an outer gap. Pairs could certainly be formed at 
the polar cap of an • B > 0 Crab pulsar but it would 
be radio-quiet or have emission very different from that ob¬ 
served. 

The growth-rate exponent described in Section 3.2 can 
also be evaluated for MSP under the assumption that Tpc = 
Ts- The distributions of A for the two sets of MSP are each 
quite compact: for the Fermi LAT set the mean value is, 
— 0.28 and for the Kramer et al set, 0.24. These 
low values in radio-loud MSP might be thought a problem, 
but as noted in Section 3.1, the magnetic flux density in the 
expression A oc is the real (unknown) value, not 

Bd- Also, MSP almost certainly have a history different from 
the typical isolated pulsar, making the implicit assumption 
of a common doubtful. 

Kramer et al find that their MSP observations are con¬ 
sistent with a single emission altitude and a compact source, 
equivalent to a light travel distance of no more than 2.4 x 10® 
cm. Pairs are not made at the polar cap of an MSP for the 
reasons given in Section 2.2 and in (i) of the Section 3.1. 
The plasma in the $1 • B < 0 case has a negligible electron- 
positron component; there is no mode suppression and there¬ 
fore nulls should not be a feature of MSP. 

But there is no reason to suppose that MSP with fl ■ 
B > 0 should not exist. (The Fermi LAT catalogue contains 
a number of non-MSP sources that are radio-quiet.) The 
problem is why are radio-quiet MSP not seen? In view of the 
estimate made by Romani (2012), this remains unresolved. 

The MSP magnetospheric structure near the light cylin¬ 
der is not well-established but we can comment that elec¬ 
trons must pass through the light cylinder in the • B < 0 
case with an integrated flux of the order of the Goldreich- 
Julian magnitude. The reason is that there exists the return- 
current problem: the star and magnetosphere inside the light 
cylinder must maintain neutrality or a constant net charge 
over time intervals of the order of the rotation period. Al¬ 


though this paper is not directly concerned with 7 -emission, 
an outer gap sustained by photon-photon collisions should 
be possible for either sign of • B. 


3.4 Nulls mode-changes and RRAT 


This class of activity is clearly a feature of increasing polar 
magnetic flux density, as demonstrated in Fig.l. In the case 
of nulls, observable radio emission ceases abruptly, within 
a time interval shorter than the rotation period, although 
in some cases very weak but detectable emission remains 
(Esamdin et al 2005, Wang et al 2007): mode-changes are 
similarly prompt. 

Pulsars listed in Tables 1 and 2 of Wang et al consist 
of 63 that have a finite null fraction, of which 5 also exhibit 
mode changes. There are also 2 showing mode changes but 
with no measured null fraction. The mean value of B is sig¬ 
nificantly larger than for the whole population and this is ap¬ 
parent in Fig. 1 in which those with mode changes are shown 
separately. The age distribution extends typically over the 
interval 4 — 40 Myr and there is some evidence (see Wang 
et al) that the null faction increases with age. 

Mode changes are much less well catalogued, but in a 
broad sense, nulls and mode-changes appear as different as¬ 
pects of the same phenomenon. A specific pulsar which ex¬ 
emplifies this is B0943-I-10 (Hermsen et al 2013, Mereghetti 
et al 2013, Bilous et al 2014). The mechanism underlying 
this phenomenon follows directly from Sections 2.2 and Sec¬ 
tion 3.1. We suggest that it is simply the suppression of 
one (or more) Langmuir modes by a change in the density 
or momentum spectrum of the electron-positron component 
on flux lines emerging from some section of the polar-cap 
area, whilst other modes remain and are observable. There 
are two obvious characteristic times which can be significant: 
the first is Tp 1 s which is the basic time-scale for fluctu¬ 
ations in electric field and plasma composition at the polar 
cap. The second is r^i, the time in which ions equivalent 
to one radiation length are accelerated from the polar-cap 
atmosphere. 


Trl KS 2.1 X 10'* 


—Psecijj 


ZsB^i\r,{\2ZV'^B-:^l‘^'' 


( 6 ) 


where ip is the angle between Cl and B. This represents the 
medium time-scale of fluctuations in Zs (see Jones 2013a). 
It is a consequence of the fact that the reduction in nuclear 
charge from Z^ to Zs by giant-dipole state formation and 
decay occurs at a limited interval of depth near the shower 
maximum. The change in Zs is unlikely to be as rapid as 
the observed onset of a null, but the change in the exponent 
A given by equation (3) need not be large to suppress the 
amplitude of the mode concerned. 

Pair production by photon-photon scattering is not es¬ 
sentially B-dependent, but the threshold momentum for an 
incoming curvature-radiation photon is fcc = 2rrJPl{khb + 
kf,b\\) and is most unlikely to be reached. There is no ev¬ 
idence that values of either Tpc or Ts are large in pulsars 
exhibiting nulls or mode changes. 

Nulls of very short duration, usually a single period, are 
frequent but we believe them to be no more than a direct 
consequence of the fluctuations with time-constant Tp that 
are a feature of the model polar cap (Jones 2013a). Times 
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given by equation ( 6 ) are not inconsistent with what are 
termed medium time-scale nulls and with mode changes. A 
null fraction increasing with age is clearly related to the age- 
related decrease in A. The cut-off value Ac is not necessarily 
well-defined because it must depend on other factors which 
have been neglected in the scaling expression for A. An ex¬ 
ample is Zs which can change in a particular region of the 
polar cap as a consequence of the medium-term instability 
with time-scale th- 

Our final comments in this Section concern the RRAT. 
These have an age distribution which differs little from the 
Wang et al nulls. But the magnetic fields are an order of 
magnitude larger. They are viewed here as merely an ex¬ 
treme case of the nulling phenomenon. They are radio-quiet 
except for intervals of one or possibly two pulses, which are 
of the order of Tp in length. We believe it is a natural con¬ 
sequence of the increased field that, for the reasons we have 
already described previously in Section 3.1, the Langmuir 
modes are suppressed except for short time intervals of the 
order of Tp. They appear to have no observable successor 
population in the ATNF logarithmic P — P diagram and, as 
we would expect on the basis of Section 3.1, are presumably 
evolving towards a radio-quiet state. 


3.5 HBRPP, CCO, XINS, AXP and SGR 

These groups of neutron stars are presently observed to be 
radio-quiet except for the HBRPP and the AXP J1810-197. 
The latter source does have coherent emission (Serylak et al 
2009) but not with the usual spectral properties, and falls 
into the class we have defined in Section 1 as secondary emis¬ 
sion. The HBRPP are defined as those with Bd > 10^® G. 
All groups except the Central Compact Objects (CCO) are 
positioned on the right-hand side of Fig. 1. Surface temper¬ 
atures are typically ~ 100 eV for the XINS, 100 — 200 eV 
for HBRPP (Kaplan & van Kerkwijk 2009). The AXP and 
SGR have much higher values, typically ~ 300 eV. 

An immediate question must be why are HBRPP radio- 
loud but the XINS, with similar values of Bd and Ts, radio¬ 
quiet? Some caution is necessary here because the XINS 
are small in number and have been detected by their 47 r X- 
ray emission: narrow radio-frequency beams may have been 
missed. It is also possible that they are all • B > 0 neu¬ 
tron stars. But for present purposes we shall discount these 
possibilities. 

Reference back to Section 3.1 indicates that coherent 
emission has two requirements. The mode growth-rate given 
by equation (3) must be sufficiently large to produce non¬ 
linearity and turbulence, and mode suppression through pair 
creation at large B must be absent. The growth-rate ex¬ 
ponent reduces to A oc for the approximation 

Tpc = Ts, following the discussion given in Section 3.2. 
Typical HBRPP surface temperature are in the interval 
100 — 200 eV and justify this approximation. The distri¬ 
butions of B for the two groups are very similar. Then for 
the XINS and for the HBRPP listed as radio-loud by Kaplan 
& van Kerkwijk, the values of A fall into two well-separated 
groups. The average values are {B]/^P~^^^) hbrpp = 7.4 
and {B]/^P~^^'^)xins ~ 1.6. This is a large difference in 
an amplitude exponent for which we have assumed that a 
value A = 30 gives adequate gain, and appears a plausible 


explanation for the presence or absence of coherent emission. 
Our only reservation here is that a common value of is 
impicit in this growth-rate comparison. The two groups may 
have different histories, with different surface properties as 
a consequence. 

We have also to comment on the CCO, in particular, the 
group of J0821-4300, J1210-5226, and J1852-b0040 observed 
by Gotthelf et al (2013). Here, the assumption Tpc = Ts 
appears satisfactory, and the values of are, 0.51, 

0.48, and 0.54 respectively. These appear to lie just above 
the notional cut-off curve that we referred to in section 3.2 
but, as in the case of the MSP, the expression requires B and 
not Bd- Thus the lack of sufficient mode growth may explain 
the absence of coherent emission now and in the future, but 
it remains possible that the CCO are f2 • B > 0 neutron 
stars. 

The magnetars (AXP and SGR) have much larger val¬ 
ues of B and of Ts- Pair creation in the special case of 
B > Be has been described by Thompson (2008). But in 
these groups, the value of A given by equation (3) is imma¬ 
terial. Following the arguments given in Section 3.1 (ii), the 
pair creation rates are much larger than in any other group. 
Mode growth is not possible. The magnetars have no obvi¬ 
ous radio-loud successors present in the ATNF logarithmic 
P — P diagram and we have to presume that, with cooling, 
mode suppression remains until A is too small to produce 
coherent emission. 


3.6 Neutron stars with $7 • B > 0 and unresolved 
problems 

There appears to be no reason why neutron stars with 
n • B > 0 should not exist. Many authors have assumed 
them to be the canonical pulsar. But it has proved diffi¬ 
cult to relate the properties of a secondary electron-positron 
plasma to the phenomena that are observed. Specifically, the 
spectral properties of the emission are essentially unchanged 
over three orders of magnitude in P and five orders of mag¬ 
nitude in Bd, whereas pair creation in many pulsars through 
conversion of curvature radiation appears impossible unless 
an ad hoc flux-line curvature is assumed (see Hibschman & 
Arons 2001, Harding & Muslimov 2002). In some recent LO- 
FAR observations (Hassall et al 2012) emission heights and 
source size do not correspond with what would be expected 
for an electron-positron plasma (Jones 2013b). 

The view of this paper is that the primary emission 
arises from the more complex plasma of 17 • B < 0 neutron 
stars. In these, the polar magnetic flux density has no essen¬ 
tial part in the emission process except for the generation 
of the acceleration field and in suppression of the Langmuir 
mode at large B. Neutron stars with 17 • B > 0 may be 
sources of secondary emission, possibly observable through 
the SKA, but otherwise emit only in the incoherent regions 
of the spectrum. Young neutron stars with 17 ■ B > 0 are 
sources of X-ray and 7 -ray emission from the vicinity of the 
light cylinder but are expected to be radio-quiet and to re¬ 
main so. A considerable fraction of Fermi LAT pulsars are in 
this category. Of 77 pulsars not classified as MSP, Abdo et 
al list 35 as radio-quiet. The radio-quiet neutron-star groups 
in this Section could have either sign of 17 ■ B. 

Two obvious unresolved problems remain within Sec- 
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tions 3.3 and 3.4. We are unable to comment here on the 
absence of radio-quiet MSP in the Fermi LAT catalogue, 
but there is some obligation to mention the well-established 
properties of long-term nulls in three pulsars; B1931-I-24 
(Kramer et al 2006), J1832-I-0029 (Lorimer et al 2012) and 
J1841-0500 (Camilo et al 2012). Whereas the time constants 
Tp and Tri have values consistent with possible relationships 
with short and medium-term nulls, times of the order of 
10 ^ — 10® s appear unrelated to any physical quantity, except 
possibly times derived from the transverse diffusion coeffi¬ 
cient of protons or ions in the polar-cap LTE atmosphere. 
A further property of these pulsars is that their spin-down 
torque is approximately halved during their null intervals. 
Although measurements of this kind are not possible in most 
pulsars, it is known that large changes are not universal. In 
one case (B0823-I-26; Young et al 2012) the fractional change 
in spin-down torque has an upper limit of 0.06. 

The likely explanation for this behaviour is that the flux 
and nature of particles crossing the light cylinder change 
during the null. The transition between ion-proton and 
electron-positron plasma is an obvious example. Our argu¬ 
ment in Section 3.2 that a state of pair production in an 
• B < 0 neutron star at Y > Xc would always be unstable 
against transition to an ion-proton plasma might be incor¬ 
rect. If this were so, such transitions could be significant be¬ 
cause all three pulsars have X ~ Xc- It is by no means clear 
which state would produce the larger spin-down torque. 


4 CONCLUSIONS 

The ideas summarized in Section 2.1 are based on well- 
established nuclear physics and in our view there can be 
no serious doubt that an ion-proton plasma is produced at 
the polar cap of an O ■ B < 0 neutron star. The motivation 
for the series of papers cited there was the belief that the 
problem of plasma composition must be solved first. Section 
2.2 summarizes ideas based on this which might be more 
contentious. The development of turbulence in the inhomo¬ 
geneous system described there is a formidable problem. The 
only simple feature is the 5-function form of velocity distri¬ 
bution which at once suggests Langmuir modes. The ques¬ 
tion of whether prior damping of the Langmuir mode might 
preclude the development of turbulence has not been ad¬ 
dressed in detail, although there appears to be no obvious 
mechanism. 

The necessary involvement of condensed matter leads 
to complexity and, in particular, to a dependence on pa¬ 
rameters that are now, and may remain in the immediate 
future, not well-known. Almost nothing is known directly 
about the surface atomic number and there is the dif¬ 
ficulty that our expression for the growth-rate exponent A 
requires the actual B which may differ from Bd- Most MSP 
are in binaries, with a history very different from isolated 
neutron stars. Nothing is known of the history of the XINS. 
The surface temperatures Ts in normal pulsars are usually 
too low to measure, the black-body photons being in an en¬ 
ergy interval strongly absorbed by the interstellar medium. 
With the complexity of the model, these factors restrict the 
nature of the predictions that can be made. 

The aim of this and previous papers has been to pro¬ 
vide a clear physical basis for understanding the nature of 


the pulsar inner magnetosphere, but only at the broadest 
level and in a qualitative way. Given the nature of the sys¬ 
tem, it may prove difficult to do more in the foreseeable 
future. Nonetheless, within these limitations, the work has 
provided an understanding of how a plasma is created that 
has the capacity to generate coherent emission over the very 
large intervals of B and P that are observed, and how that 
emission eventually ceases. The finding that relatively small 
(background) fluxes of electrons and positrons can affect 
mode growth has proved very significant. It is likely to be 
the basis for the magnetospheric changes that are seen as 
nulls and mode changes. It also demonstrates why some very 
active neutron stars have no observable coherent emission. 
More observational evidence on source sizes and altitudes 
would be of interest, as would high-sensitivity SKA stud¬ 
ies of what we have referred to as secondary coherent emis¬ 
sion. An understanding of the inner magnetosphere region of 
high energy densities may be an essential first step towards 
achieving a more complete understanding of the whole mag¬ 
netosphere. 
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